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Abstract 

The phase behavior of asymmetrically interacting ABC star triblock copolymer melts is in¬ 
vestigated by the self-consistent field theory (SCFT). Motivated by the experimental systems, 
in this study, we focus on the systems in which the Flory-Huggins interaction parameters sat¬ 
isfy Xac > Xbc ~ Xab- Using various initialization strategies, a large number of periodic struc¬ 
tures have been obtained in our calculations. A fourth-order pseudospectral algorithm com¬ 
bined with Anderson mixing method is used to compute the free energy of candidate structures 
carefully. The stability has been detailedly analyzed by splitting the free energy into internal 
and entropic parts. A complete and complex triangular phase diagram is presented for a model 
with Xac > Xbc = Xab in which fifteen ordered phases, including two-, and three-dimensional 
structures, have been predicted to be stable from the SCFT calculations. Generally speaking, 
with the asymmetrical interactions, the hierarchical structures tend to be formed near the B- 
rich corner of the triangular phase diagram. This work broadens the previous theoretical results 
from equal interaction systems to unequal interaction systems. The predicted phase behavior 
is in good agreement with experimental observations and previous theoretical results. 

*To whom correspondence should be addressed 
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Introduction 


Block copolymers, constructed by linking together ehemieally distinet subehains or bloeks, spon¬ 
taneously assemble into exquisitely ordered soft materials .^The self-assembled order struetures, 
spanning length seales from a few nanometers to several mierometers, offer a diverse and expand¬ 
ing range of praetieal applications in, for example, optical materials, mieroeleetronie materials, 
drug delivery, advanced plasties, and nanotemplates 

The development of nanotechnology using block copolymers requires a good understanding of 
the phase behavior of the bloek eopolymers. The self-assembling meehanism of bloek eopolymers 
sensitively depends on bloek types, the number of bloeks, bloek-bloek interaetions, arehiteeture, 
and topology of the block polymers. The equilibrium ordered patterns ean be formed due to the 
delicate balanee between these eompeting faetors. In ABC triblock copolymers, the number of 
controlled parameters is at least five, including three interaction parameters XabN, XacN, XbcN, 
and two independent bloek eompositions /a and /g. A is the degree of polymerization and Xap 
is the Flory-Huggins interaetion parameter eharaeterizing the interaetion between two ohemieally 
different bloeks a and /3. Compared with linear copolymers, the star-shaped eopolymers have 
complicated phase behavior or physieal properties indueed by different moleeular arehiteeture. In 
the ordered phases of ABC star tribloek eopolymers, the most distinet feature is the arrangement of 
the junetion points. If the ehain lengths of three bloeks are eomparable, junetion points are aligned 
on a one-dimensional (ID) straight line, then eylindrieal morphologies ean be formed naturally. 
Their eross seetions tend to show two-dimensional (2D) patterns sinee polymer/polymer interfaees 
can be flat surfaees due to the repulsion forces between “unconneeted” branehes. These faetors 
lead to the formation of 2D polygonal tiling patterns, or Arehimedean tiling. The tiling patterns 
ean be eneoded by a set of integers [k, /, m,... ], indieating that a k-gon, an l-gon, and an m-gon, 
ete., meet eonseeutively at eaeh vertex. If some asymmetry, to the eontrary, is introdueed to the 
eompositions, the junetion points are mostly aligned on the eurved trails. Consequently three- 
dimensional (3D) structures ean be formed. Furthermore, when the interactions are strong enough, 
the ABC star triblock copolymers can self-assemble into hierarchical structures 
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Owing to the rich phase behavior, a number of experiments have been carried out on the 
phase and phase transition of ABC star triblock copolymers in past decades .ISHIH particular, 
Matsushita and co-workers have conducted systematic studies on the morphologies formation 
of (polyisoprene-polystyrene-poly(2-vinylpyridine)) (ISP) star triblock copolymers. Several or¬ 
dered tiling patterns, such as [6.6.6], [8.8.4], [12.6.4], [3.3.4.3.4], and even a dodecagonal sym¬ 
metric quasicrystalline tiling ,2^ have been observed in the ISP star triblock copolymers .ISEEUMI 
Meanwhile, several hierarchical structures including cylinders-in-lamella, lamellae-in-cylinder, 
lamellae-in-sphere, and hierarchical double gyroid structures are also discovered with the asym¬ 
metries of composition .IhEEl Recently, Nunns et al. used polyisoprene (I), polystyrene (S) and 
poly-(ferrocenylethylmethylsilane) (F) to synthesize ISF star triblock copolymers ^ Three 2D 
patterns, including [8.8.4], [12.6.4], and lamellae with alternating cylinders, were observed. Both 
of the experimental systems satisfy asymmetrical interaction between different blocks, i.e., XacN > 
XbcN ~ XabN. In particular, in ISP star triblock copolymers, the interaction strengths are known 
to follow the order Xis^ ~ Xsp^ < Xip^ while in ISF star-shaped triblock copolymers, XsF^ ~ 
XisN < XifN 

Besides experimental works, theoretical studies provide a good understanding of phase be¬ 
havior of ABC star triblock copolymers. Bohbot-Raviv and Wangle used a coarse-grained free 
energy functional to numerically investigate some morphologies of ABC star triblock copolymers. 
In 2002, Gemma and co-workers carried out Monte Carlo (MC) simulations on ABC star tri¬ 
block copolymers with equal interactions between the three components. The phase behavior of 
ABC star triblock copolymers with composition ratio of, /a : /b : /c = TO : 1.0 : x, was investi¬ 
gated in detail in strong segregation region. Five kinds of 2D cylindrical phases, three kinds of 
lamellar-type phases and two kinds of continuous matrix phases were obtained from the MC sim¬ 
ulations. Huang and co-workers studied the effects of composition and interaction parameter on 
the phase behavior of ABC star copolymers with equal interactions among the three components 
using dissipative particle dynamics (DPD) simulations. Several efforts have also been made to 
study the phase behavior of ABC star triblock copolymers using the SCFT .ISHUI 2004, Tang 
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et al. ^ started to use a 2D SCFT simulation to study the phase behavior of ABC star triblock 
copolymers. Based on the SCFT, Zhang et al. 123 and Li et al. 123 examined the weak and interme¬ 
diate segregation cases of mainly 2D structures with equal interactions. Zhang et al. l23also chose 
XabN = XbcN = 25.0, XacN = 37.0 to model the ISP star triblock copolymer system of the type 
Ai.oBi.oCx. Accordingly, a ID phase diagram was obtained as a function of x from 0.5 to 2.0. The 
stability of the different lamellar morphologies formed from ABC star triblock copolymers has 
been examined by Xu et al. ^ These SCFT studies mainly focus on the symmetrically interacting 
systems and the phase behavior of 2D structures. Despite these previous experimental and theo¬ 
retical studies, a comprehensive understanding of ABC star triblock copolymers is still lacking, 
especially for the asymmetric interaction systems. In this work, we will explore the phase behav¬ 
ior of ABC star triblock copolymers, including 2D and 3D structures, with asymmetric interaction 
parameters between chemically different blocks. 

Theoretical approaches to investigating the phase behavior of block copolymers, often involve 
minimizing an appropriate free energy functional of the system, and comparing the free energies of 
different candidate structures. Therefore a systematic examination of the emergence and stability 
of ordered phases requires the availability of suitable free energy functionals and accurate methods 
to compute the free energy of ordered phases. Owing to a large number of studies, it has been well 
proven that the SCFT provides a powerful theoretical framework for the study of the phase behavior 
of block copolymers .ESEl! jn particular, the SCFT can be used to determine the relative stability 
of different phases because it provides an accurate estimate of the free energy. The essence of the 
SCFT is that the free energy of the system can be written as a functional of the spatially varying 
polymer densities and a set of conjugate fields. Minimizing the free energy functional with respect 
to the densities and conjugate fields leads to a set of equations, encoded as SCFT equations. The 
SCFT equations are a set of highly nonlinear equations with multi-solutions which correspond to 
the different ordered phases of block copolymers. The equations also have a strong nonlocality 
that emerges from the connection of propagators and densities, conjugate fields. 

Solving the SCFT equations requires iterative techniques. Owing to the nature of iterative 
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methods, the solutions sensitively depend on the initial eonfiguration at the start of the iteration. 
A series of effieient strategies of sereening initial eonditions are developed based on the faet that 
all periodie struetures belong to one of 230 spaee groups p^ l 30|3i | iterative step, effieient 

numerieal sehemes are required to solve the propagator equations. In the past years, two eomple- 
mentary methods, ineluding speetral methodsand real-spaee methods!^ to solve the SCFT 
equations have been developed. In reeent years, an effieient pseudospeetral method has been intro- 
dueed to solve propagator equations .E2H1I! algorithm takes advantage of the best features of 
real- and Fourier-spaee with the eomputational effort seale of 0{NlogN) based on the fast Fourier 
transformation (FFT). N is the number of speetral modes or diserete points in real-spaee. To obtain 
the solutions of the SCFT equations eorresponding to the saddle points of the SCFT free energy 
funetional, the iterative methods are required to make the iteration eonvergent. The quasi-Newton 
methods were employed in the fully speetral approaeh to the SCFT by Matsen and Sehiek .1^ A 
simple mixing method by a linear eombination of two eonseeutive fields was introdueed by Drolet 
and Fredriekson^for the SCFT simulations. Reeently the Anderson mixing methodhas been 
proven by itself to greatly reduee the number of SCFT iterations. From the perspeetive of nonlinear 
optimization, Cenieeros and Fredriekson^ devised a elass of effieient semi-implieit sehemes for 
solving the SCFT equations using the asymptotie expansion teehnology. Later, Jiang et al. I^have 
extended these algorithms to the SCFT ealeulations for multieomponent polymer systems. 

A generie strategy of theoretieal studying phase behavior of eomplex bloek eopolymer systems 
ineludes two steps .l ^^ l ^^ lThe first step involves an effieient strategy to produee a library of possible 
eandidate struetures. In the seeond step, the eandidate struetures are used as initial eonditions in the 
more aeeurate methods to eompute free energies, whieh are used to eonstruet phase diagrams of the 
systems. In this work, we apply this strategy to examine the phase behavior of ABC star tribloek 
eopolymers using the SCFT. Speeifieally, the strategies developed in our previous work are used 
as a sereening teehnique to obtain eandidate struetures as many as possible. These strategies in- 
elude (1) knowledge from previous experiments and theories; (2) knowledge from related systems, 
for example, dibloek eopolymers; (3) eombination and interpolation of known struetures; and (4) 
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random initial configurations. Using these candidate structures as initial conditions, a fourth-order 
pseudospectral method combined with Anderson mixing method is employed to study the stability 
of ordered phases. To model the asymmetric interacting experimental systems of ISP and ISF, the 
interaction parameters of XacN > XbcN = XabN are used in our study. It should be emphasized 
that, to broaden the scope of the research, the 2D and 3D ordered phases are included in our cal¬ 
culations. The resulting free energies of the different ordered phases are used to construct phase 
diagrams. 

Theory and Methods 

We consider an incompressible melt of n ABC star triblock copolymers with a degree of poly¬ 
merization A in a finite volume of V. The chain lengths, or compositions, of A, B, and C blocks 
are /a A, /gA, and fcN (/a +/s + /c = 1)> respectively. A characteristic length of the copolymer 
chain can be defined by the radius of gyration, which is used as the unit of length, so that all spatial 
lengths are presented in units of Rg. Within the mean-field approximation to statistical mechanics 
of the Edwards model of polymers, at a temperature T, the free energy functional F per chain of 
the Gaussian triblock copolymer melt is 

-^^-logQ+^ fdri Xap^^a(r)^p(r)-J^Wa(r)Mj^)-ri(r)ll-J^^a(r)]}, (1) 

riKBi V J i- Q, a ^ 

where a, /3 G {A, B, C} are the block labels, is the monomer density of the a-blocks, and Q is the 
partition function of one star block copolymer chain in the mean field, Wo:(r), produced by the sur¬ 
rounding chains, rj (r) is the Lagrange field to ensure the local incompressibility. The interactions 
between the three chemically distinct monomers are characterized by three Flory-Huggins inter¬ 
action parameters multiplied by polymerization degree, i.e., XbcN, XacN, and XabN. First-order 
variations of the free-energy functional with respect to the monomer densities and the conjugate 
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fields subjected to the incompressible condition lead to the following set of SCFT equations 


WA (r) = Zab^^b (r) + ZAcN^cir) +1? (r), 

( 2 ) 

vvs(r) = XABN^Air) + XBcN^c{r) +1? (r), 

(3) 

wc(r) = XacN 4 >a (r) -f XbcNM^) + ^ (r), 

(4) 

1 = ())A(r)-f <))s(r) + <))c(r), 

( 5 ) 

1 rfA , 

^Air) = - dsqA{r,s)q\{rjA-s), 

Q ./o 

( 6 ) 

1 ffs 

= qI^ dsqB(r,s)ql^(rjB-s), 

( 7 ) 

1 ffc 

<^>c(r) = gy^ dsqc{r,s)q]^{rJc-s), 

( 8 ) 

Q= ^ J drqK(r,s)ql(rjK-s), ^s^K. 

( 9 ) 


In these expressions, the functions qK{r,s) and <5r^(r,5) {K & {A,5,C}) are the end-integrated seg¬ 
ment distribution functions, or propagators, representing the probability of finding the 5-th segment 
at a spatial position r. These propagators satisfy the modified diffusion equations (MDEs) 


j-^qK{r,s) = VlqK{Y,s)-WKqK{r,s), ^x(r,0) = 1, 5 G [OJk], 
Y^^l{r,s) = Vlql{r,s)-WKqlir,s), ^l(r,0) = qL{rJi)qM{rjM), 


( 10 ) 


where (KLM) G {(ABC), (BCA), (CAB)} . Numerically solving these SCFT equations involves 
an iterative procedure starting with an initial configuration of the fields wk(y) with K G A^B^C. 
The MDEs (Eqn. (??)) are then solved to obtain the propagators, which are used to compute the 
densities ^k{^) and to update the mean fields wk(y). The iteration is continued until these mean 
fields and densities are self-consistent such that they satisfy the SCFT equations within a prescribed 
numerical accuracy. 

There are two main steps in studying the phase behavior of block copolymers. The first step 
is to obtain possible candidate structures as many as possible. Several strategies of exploring or- 
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dered phases have been proposed in our previous works Beyond the random initial values, 
these approaehes include (1) knowledge from experiments and theories, such as small-angle X-ray 
scattering images, space group theory for periodic structures; (2) knowledge from related systems; 
(3) combination and interpolation of known structures. Using these diverse strategies of initializa¬ 
tion, a large number of ordered phases can be generated as the solutions of the SCFT equations. 
In addition, the possible candidate structures from previous theoretical and experimental studies 
are considered in our calculations. The second step is to identify the stability of these phases by 
comparing their free energies using efficient numerical methods. In this work, we combine an 
improved pseudosepectral method with Anderson mixing algorithm to solve the SCFT equations 
for periodic block-copolymer morphologies. A fourth-order accurate Adams-Bashford scheme^ 
is used to discretize the MDEs. The initial values required to apply this formula are obtained 
using a special extrapolation method based on the second-order operator-splitting scheme 
A modified integral formula for closed interval is chosen to solve the integrated equations (??)- 
(??) that can guarantee fourth-order precision in 5-direction whether the number of discretization 
points is even or odd To ensure the accuracy, we require that these substeps of contour length 
s are smaller than 0.01 in the fourth-order accurate scheme. The FFT is used to translate the data 
between real- and Fourier-space in the pseudospectral method. 

To obtain the equilibrium morphologies of SCFT equations, iterative methods shall be required 
to update the conjugate fields. For this step we choose the Anderson mixing algorithm, firstly 
proposed by Anderson ,1^ then introduced into polymer theoretical calculations by Schmid and 
Muller ,1^ Thompson et al. Owing to the local convergence of this method, we use the simple 
mixing method alone at the start of the algorithm to obtain better initial values for the fields, fol¬ 
lowed by Anderson mixing approach alone to accelerate the convergent procedure to the prescribed 
accuracy in the fields. The Anderson mixing method requires the fields of previous n steps when 
update new fields. Using the previous n step fields, the Anderson mixing method produces an 
R-order linear equations system from a least square problem. Then the new fields will be obtained 
by a combination of the fields of previous n steps, with the solution of linear system being the 
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weight factors. For relatively simpler systems, such as diblock copolymers, the Anderson mix¬ 
ing algorithm can significantly reduce the required number of iterations with few histories .EUSI 
For more complex situation, a larger n shall be taken to update fields to accelerate the convergent 
procedure. We find that assembling the n-order linear system spends more computation time than 
solving this linear system. When n becomes too large, it will slow the iteration. Here we overcome 
this problem by rearranging the elements of the nxn matrix. The technical details can be found in 
the Appendix section. By using our approach, the Anderson mixing method can robustly converge 
without slowing the SCFT iterations. In practice, we use the available histories of 50 steps. 

Here we only consider the periodic structures, therefore, periodic boundary conditions are im¬ 
posed on each direction. In our calculations, all spatial functions are all expanded in terms of plane 
waves. For 2D morphologies, a square box is simulated. 128 x 128 plane-wave basis functions are 
used to discretize the 2D box. For 3D structures, we use a cubic unit cell in most of our calcu¬ 
lations. The number of plane-wave basis functions is 32 x 32 x 32. The size of computation box 
plays an important role in determining the stability of ordered phases. For a given phase, its free 
energy is minimized with respect to the box sizes by the steepest descent approach coupled with 
solving the SCFT equations 

Based on the discussion above, we summarize the iteration procedure by sketching the numer¬ 
ical recipes: 

Step 1 Starting from the given initial conditions and computational box, the fourth-order pseu- 
dosepectral method combined with Anderson mixing method is applied to obtain the ordered 
phases when the field’s change is smaller than the prescribed error Ci. 

Step 2 Optimizing the size of unit cell by minimizing the free energy with the steepest descent 
method. 

Step 3 Goto Step 1 until the free energy change is smaller than a given error £ 2 . 

To ensure enough accuracy, in our implementation, each calculation is terminated until the field’s 
change (defined in Appendix) at each iteration is reduced to £1 = 10^^ (corresponding to a free 
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energy ehange of about 10^^), and £2 = 10^^. For the eases of XacN = 50.0, XabN = XbcN = 30.0 
in the eurrent work, the relative values of free energy among different phases in determining the 
phase boundary are from 0(10^^) to 0(10^^). Therefore our numerieal resolution is adequately 
aeeurate for eonstrueting phase diagrams. 

Results and Discussion 

The previously theoretieal studies mainly foeus on the equally interaeting ABC star tribloek eopoly- 
mer systems, i.e., Xbc^ = XacN = XabN ^ | 22 | 24 | 25 ] 27 | pjo^rever, under experimental eireumstanees 
as mentioned above, many ABC star tribloek eopolymers with asymmetrie interaetions have been 
synthesized to observe their phase behaviors, sueh as ISF^ and ISpl^l star tribloek eopolymers. 
In partieular, the interaetion parameters of the above systems satisfy the relationship XabN ~ 
XbcN < XacN. In order to make a meaningful eomparison of the theoretieal study and exper¬ 
imental investigation, the interaetion parameters are ehosen in the ealeulations so that they are 
appropriate for these experimental systems. Although aeeurate values of the Flory-Huggins pa¬ 
rameters are not available in the literatures, qualitative interaetion strengths are known to follow 
the order Xis ~ Xsp < Xip ^^Xsf ~ Xis < Xif -^In what follows we ehoose asymmetrie interae¬ 
tion parameters, XacN = 50.0, XabN = XbcN = 30.0, and equal statistieal segment lengths. This 
is a rough approximation to ISP and ISF eopolymers, with the important differenee that, in the 
idealized system, we negleet the small differenees between XabN and XbcN. 

As mentioned above, there are two main steps in studying the phase behavior of bloek eopoly- 
mers. The first step is to obtain as many possible eandidate struetures as possible. The seeond step 
is to identify the stability of these patterns by eomparing their free energies and eonstruet the phase 
diagrams. In order to further analyse the stability of eandidate patterns, it is helpful to split the free 
energy into two parts: internal {U) and entropie {—TS). The internal and entropie eontributions to 
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the free energy ean be expressed as 


■»£r = 


( 11 ) 


where a,/3 G {A,B,C} are the bloek labels. 


Candidate Patterns 

Using different initialization proeedures, a large number of eandidate struetures have been ob¬ 
tained in our previous work .I^Here we only present the stable phases in the ease of XacN = 50.0, 
XabN = XbcN = 30.0, as shown in Fig.[^ Among these eandidate phases, several polygon tiling 
patterns, or the eylindrieal struetures, with translation invariant along the third direetion are found 
as 2D phases. These polygon patterns inelude [6.6.6], [8.8.4], [12.6.4], [3.3.4.3.4], [8.6.6;8.6.4], 
[10.6.6;10.6.4;8.6.6;8.6.4]. The first pattern is designated as [6.6.6] beeause eaeh vertex in this 
tiling is surrounded by three hexagonal polygons. Similarly the seeond and third patterns are 
named [8.8.4] and [12.6.4], respeetively. The fourth pattern looks more eomplex. There are two 
types of vertiees: one is surrounded by 10-gon, 8-gon and 4-gon, whereas the other is formed 
by a deeagon, a hexagon and a tetragon. From our naming eonvention, it should be named 
[10.8.4;10.6.4]. However, with the help of triangles and squares of hypothetieal tiling, as Fig.[2 
shows the superimposed tiling on a sehematie drawing, it is noted all the meeting vertiees are sur¬ 
rounded by three regular triangles and two squares, whieh is one of the Arehimedean tilings. In 
order to eompare with experimental results, we eneode it as [3.3.4.3.4] .I^The fifth pattern also 
possesses two kinds of vertiees: one is surrounded by 8-gon, 6-gon, and 4-gon, whereas the other is 
formed by an 8-gon and two 6-gons. Consequently, this pattern is designated as [8.6.4;8.6.6]. Sim¬ 
ilarly, the fifth pattern has four ki nds of vertiees, therefore, it is named [10.6.6;10.6.4;8.6.6;8.6.4]. 
Besides these polygonal phases, additional three 2D struetures, i.e., three eolor lamellae (LAM), 
the eore-shell eylinders (HC), and hierarehieal eylinders-in-lamella phases (L-i-C) are also obtained 
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Figure 1: Ordered phases of ABC star triblock copolymers obtained using the SCFT calculations 
with XacN = 50.0, XabN = XbcN = 30.0. The colors of red, green, and blue, indicate the regions 
where the most components are A, B, and C, respectively. 


12 











in our simulations. At the same time, a series of 3D struetures are obtained. From the morphol¬ 
ogy of patterns, these 3D patterns ean be elassified into eore-shell phases and hierarehieal strue¬ 
tures. The former ineludes eore-shell spheres in body-eentered-eubie lattiee (BCC) and eore-shell 
double-gyroid phases (DG). The hierarehieal patterns eonsist of two kinds of hierarehieal eylin- 
ders paeked hexagonally (HHC), two kinds of eylinders-in-lamella phases, with the eylinders being 
paeked hexagonally (HPL) and tetragonally (TPL), and hierarehieal gyroid phases (HDG). Note 
that due to the equal interaetion parameters of XabN = XbcN, these struetures have their mirror 
phases along the phase path of isopleth /a = fc- 

Phase Behavior of A] oBi .oCx Star Triblock Copolymers 

In a series of experiments h2 | i3 | i7 | 46 ] three arms are kept to be of equal length and the 

arm-length ratio is expressed as 1.0 : 1.0 : x. Motivated by these experiments, we start with the 
ealeulation of the phase stability along this phase path. Here we assume that A and B arms have 
equal length and the C arm holds the arm-length ratio of x = fc/f a- The free energy differenees of 
eandidate phases from the value of the homogeneous phase as a funetion of the volume fraetion of 
fc are given in Fig.|^(a). The phase stability regions as a funetion of x are presented in Fig.[^(b). 
In Fig.[^(a), the free energies of some metastable phases along this path are not shown, sueh as that 
of [10.6.6;10.6.4;8.6.6;8.6.4] in the region of 0.32 < /c < 0.46, where it has higher free energy 
than [6.6.6], or [8.8.4]^, or [3.3.4.3.4]. With the inerease of x, the phase sequenee is LAM^ 
[8.8.4] 1 ^ [8.6.6;8.6.4] ^ [6.6.6] ^ [8.8.4]^ ^ [3.3.4.3.4] ^ [12.6.4] ^ TPL LAM^. The 
eorresponding stable regions are x < 0.435 (LAM^), 0.435 < x < 0.67 ([8.8.4]^), 0.67 < x < 0.71 
([8.6.6;8.6.4]), 0.71 < x < 1.04 ([6.6.6]), 1.04 < x < 1.39 ([8.8.4]2), 1.39 < x < 1.68 (3.3.4.3.4), 
1.68 < X < 2.05 ([12.6.4]), 2.05 < x < 2.76 (TPL), and x > 2.76 (LAM^), respeetively. When 
fc is small, in one period, the lamellar strueture of LAM^ of CACB type ineludes two thiek A 
and B layers, and two thin C layers. The symbols of [8.8.4]^ and [8.8.4]^ are used to distin¬ 
guish the 4-eoordinated polygons segregated by different bloek. In the [8.8.4] ^ phase, the minority 
C-arms form the 4-eoordinated domains, and bloeks A and B alternatively form 8-ooordinated 
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Figure 2: (a) Free energy differences from the value of the homogeneous phase as a function of 
the volume fraction of C composition for ABC star triblock copolymers with symmetric A and B 
arms, (b) Phase stability regions as a function of the arm-length ratio of v = /c//a if a = /s)- Note 
that in the [8.8.4] ^ phase, the minority C blocks form the 4-coordinated domains, and blocks A and 
B alternatively form 8-coordinated microdomains. While in the [8.8.4]^ morphology, the A and C 
blocks form the 8-coordinated polygons, and B blocks form the domains with 4-coordinations. In 
one periodicity, the LAM^ phase has the CACB lamellar sequence, whereas the LAM^ structure 
has the BABC layers. 
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polygons. While the unit eell of [8.8.4]^ eontains one 4-eoordinated B-domain, 8-eoordinated and 
4-eoordinated mierodomains formed by bloeks A and C, respeetively. From the phase path, in 
tiling patterns, the eoordination number of C-domains is proportional to the volume fraetion of fc- 
With the inerement of fc from 0.18 to 0.50, the eoordinations of C-domains ehange from 4, to 6, 
to 8, then to 10 (in the [3.3.4.3.4] phase), finally to 12. Further inereasing fc, the asymmetry of A 
(B) arm and C arm results in arrangement of junetions on a eurve line. Then a 3D strueture of TPL 
appears when 0.51 < /c < 0.58. When fc > 0.58, the lamellar phase LAM^ is stable in whieh the 




Figure 3: (a) Internal energy of AU/nksT and (b) entropie energy of —S/nksT 
as a funetion of fc on the phas^ath of = /s. The morphologies of [8.8 
and LAM^ are explained in Fig. [2 


of various struetures 
.4]i, [8.8.4]2, LAM^ 


arrangement manner is BABC type in one period, ineluding one thiek C layer and three thin BAB 
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layers. 

We ean analyze the stability of these ordered phases after splitting the SCFT energy funetional 
into the internal and entropie parts (see Eq. (??)). Fig.|^ gives the internal energy At/ /nksT, sub- 
traeted by that of the homogeneous phase together with the entropie energy of —S/nks, as a 
funetion of fc- From Fig.|^(a), we ean find that FAM^ (CAGE layers) has very high internal en¬ 
ergy at small fc whieh is indueed by the penetrations of A and B arms through C domains At 
the same time, the lamellae of FAM^ is favorable from the aspeet of entropie energy, beeause the 
A and B bloeks ean get the largest entropy in the lamellar strueture when they have equal large 
lengths eompared with the C bloeks. The eombination of the two eontributions makes the FAM^ 
be the stable phase when fc <0.18. When inereasing fc, the internal energy beeomes dominant 
and the stable phase transfers from FAM^ to the eylindrieal struetures where the arm penetra¬ 
tions in dissimilar phase regions are diminished. When the stable phase transfers from [8.8.4]^ 
to [3.3.4.3.4] (also termed as [10.8.4; 10.6.4]), the entropie energy plays a dominant role on the 
stability. It is attributed to the inerement of C arm whieh tends to form large C-domains. The 
enlarged C-domains have an opportunity to meet mueh more A-domains whieh will inerease the 
internal energy due to the largest interaetion parameter XacN- On the other hand, as the arm C 
inereases, the A and B arms beeome shorter. The A and B arms ean freely streteh in their rieh 
domains whieh ean greatly reduee the entropie energy. The explanation is also available to the 
appearanee of [12.6.4], TPF as stable phases. Further inereasing fc to 0.58, the lamellar phase of 
FAM^ (BABC layers) is stable again. The reason is similar to that of the stability of FAM^ phase 
when fc is smaller than 0.18. 

Along this phase path, Matsushita and eo-workersl^ synthesized a set of Ii.oSi.qPx eopolymers 
and observed a number of ordered struetures. Four ISP star tribloek eopolymers with volume ratios 
of 1.0:1.0:0.7,1.0:1.0:1.2,1.0:1.0:1.3, and 1.0:1.0:1.9 were investigated. The eylindrieal struetures 
of Ii.oSi.oPo.77 Ii.oSi.oPi.27 Ii.oSi.oPi.37 andIi.oSi.oPi.9 are [6.6.6], [8.8.4], [3.3.4.3.4], and [12.6.4], 
respeetively. From our simulations, the resulting phase behavior is in good agreement with these 
experimental measurements when 0.7 < x < 2.0. At higher asymmetries, Takano et al. I^observed 
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an L+S phase in I 1 . 0 S 1 . 0 P 0.2 system, and an L+C phase in I 1 . 0 S 1 . 0 P 3.0 and I 1 . 0 S 1 . 0 P 4.9 eopolymers. 
However, in our ealeulations, the LAM phase dominates the regions. This might be attributed to 
the relatively weak interaetion parameters of XabN or XbcN in our simulations so that a further 
segregation between eomponents A and B (or B and C) ean not oeeur. In 2007, Hayashida et al.^ 
have diseovered a eylinders-in-lamella phase in whieh the staeking manner of the eylinders seems 
to be random in the experiments of ISP star tribloek eopolymer melts. In our ealeulations, the 
eylinders-in-lamella, TPL, is found to be stable when 0.51 < /c < 0.56, in whieh the eylinders are 
staeking tetragonally. The diserepaney is attributed to the thermodynamie fluetuations whieh may 
affeet the arrangement of eylinders under experimental eireumstanees. While within the mean-field 
level theory, the fluetuations have been negleeted. 

Along the similar phase path, Nunns et al. synthesized a set of ISF star tribloek eopolymers 
and observed the polygonal tilings [8.8.4] and [12.6.4]. Our resulting phase behavior is generally 
eonsistent with the ISF experiments. A deviation between our theoretieal results and experiments 
is that the [12.6.4] was observed in I 0 . 40 S 0 . 37 F 0.23 system. The diserepaney ean be attributed to the 
different interaetions, and different monomer sizes. 

Our eomputational results also agree with the previous theoretieal ealeulations Among 

these works, Zhang et al. I^eonsidered 2D tiling patterns and ehose the asymmetrie Flory-Huggins 
interaetion parameters of XacN = 37.0, XabN = XbcN = 25.0 to model the system of ABC star 
tribloek eopolymers. A ID phase diagram of the star tribloek eopolymers Ai oBi oCx was obtained. 
With the inerease of x, the phase transition ehanges from [8.8.4]^ to [8.8.4]^, then to [8.6.6;8.6.4] 
and finally to [12.6.4]. The eorresponding stable regions are 0.50 <x< 0.86, 0.94 < x < 1.33, 
X ~ 1.45 and 1.57 <x< 2.00, respeetively. The phase behaviors of [8.8.4]\ [8.8.4]^, and [12.6.4] 
are qualitatively eonsistent with our results. There are some diserepaneies between their results 
and our eomputer simulations. In their phase diagram, the Arehimedean tiling of [3.3.4.3.4] is 
not ineluded in their simulations, the order pattern [6.6.6] is metastable, and the stable area of 
[8.6.6;8.6.4] is different from ours. The diserepaneies ean be attributed to two aspeets. The first one 
is that more eandidate struetures are involved in our simulations. The seeond one is the differenee 
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of the Flory-Huggins interaction parameters. 


Phase Behavior of Ai qBxCi o Star Triblock Copolymers 


On the above phase path, many candidate structures, such as 2D polygonal phases [10.6.6;10.6.4; 
8.6.6;8.6.4] and 3D hierarchical structures of HHC, HDG, do not appear. To obtain the stability 
regions of the 2D and 3D phases which have been observed in experiments, we turn to another 
phase path of isopleth A = fc, i.e., by fixed equal length of arms A and C and the arm-length 
ratio of 1.0 : x : 1.0 with an increment of volume fraction fg of 0.01. The free energy differ¬ 
ence from the value of the homogeneous phase as a function of fg varying from 0.20 to 0.68 
is plotted in Fig.|^(a). The phase stability regions as a function of x are presented in Fig.|^(b). 
Along this phase path, the phase sequence with increasing fg is LAM —[8.8.4] —)■ [6.6.6] —)■ 
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Figure 4: (a) The free energy differences of various candidate phases, from the value of the ho¬ 
mogeneous phase, as a function of fg for ABC star triblock copolymers with symmetric A and C 
arms, (b) Phase stability regions as a function of the arm-length ratio of x = fg/fA (fA = fc) of 
ABC star triblock copolymers with asymmetric interactions, XacN = 50.0, XabN = XbcN = 30.0 
as a function of the arm-length ratio ofx = fg/f a- 
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[10.6.6;10.6.4;8.6.6;8.6.4] ^ [12.6.4] ^ L+C HDG -> HHC. The stable region of LAM phase 
is /s < 0 .24. At the eenter part of the phase path, the ehain lengths of three arms are close to one 
another, junction points are aligned on a straight line, and hence 2D tiling patterns can be formed. 
The stability regions of 2D cylindrical phases are 0.24 < /s < 0.33 ([8.8.4]), 0.33 < /s < 0.46 
([6.6.6]), 0.46 < /s < 0.52 ([10.6.6;10.6.4;8.6.6;8.6.4]), /s 0.53 ([12.6.4]), respectively. From 
the phase path, we can also find that, in polygonal patterns, the coordinations of B-domains are 
proportional to the volume fraction of /g. With the increase of /g from 0.24 to 0.53, the coordi¬ 
nation number of B-domains goes through a gradual change from 4 ([8.8.4]^), 6 ([6.6.6]), 8 or 10 
([10.6.6;10.6.4;8.6.6;8.6.4] tiling), to 12 ([12.6.4]). 

Further increasing /g, the A and C arms become shorter. Also, owing to the largest interaction 
parameter XacN, a further segregation between two minority blocks A and C occurs within the 
large-length-scale phase, and the system can form some hierarchical morphologies. The interesting 
hierarchical patterns include 2D L-i-C phase, and 3D patterns of HDG, HHC, and their stability 
regions are 0.54 < /g < 0.63 (L-i-C), 0.63 < /g < 0.67 (HDG), /g > 0.67 (HHC), respectively. 

After splitting the SCFT energy functional into internal energy and entropic energy as ex¬ 
pressed in Eqn. (??) (see Fig.|^, the stability of 3D hierarchical structures HHC and HDG along the 
phase path can be understood more readily. When the ratio of x = /g //4 (/a = fc) is large enough, 
the asymmetric ABC star triblock copolymers will exhibit similar phase behavior of asymmetric 
diblock copolymer. The longest B arm plays an equal role as the long block in an asymmetric 
diblock copolymer, whereas the short A arm, together with C arm, are just like the short block in 
the diblock copolymer. In asymmetric diblock copolymers, the systems tend to curve the interface 
towards the minority domain. It requires the minority block to stretch, and the cost is more than 
compensated for by relaxation of the longer blocks, which increases the internal energies Be¬ 
sides the asymmetry of blocks, the segregation power between arms A and C, XacN, is sufficiently 
large to separate the A-, and C-domains, leading to the formation of hierarchical structures. In 
particular, for the asymmetric ABC star copolymers, the blocks A and C have shorter chains than 
the short block in AB diblock when the curve interface can be formed. The A and C arms can 
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Figure 5: (a) Internal energy of AU/nksT and (b) entropie energy of —S/nksT of various struetures 
as a funetion of fg on the phase path of /a =fc- 
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freely relax in their paeking frustration domains whieh greatly reduees the entropic energy of the 
system. Although the internal energies of 3D hierarehical struetures are higher than that of L+C 
structure, their entropic energies are much lower than that of the 2D cylinders-in-lamella phase. As 
a consequence, the combination of the two competing energies makes 3D hierarchical structures 
stable rather than L+C when fg >0.63. At higher asymmetries, block copolymers prefer to form 
structures with a higher interfacial curvature Therefore HHC phase with larger spontaneous 
curvature than that of HDG phase, has lower free energy when fg > 0.67. 

Triangular Phase Diagram 

The phase transition sequence for systematically varying volume fractions can be obtained by 
repeating the free energy comparison among the candidate structures. The results of the phase 
transition sequences can be summarized in terms of phase diagrams. For the case of asymmetric 
interaction parameters of XacN > XabN = XbcN, the triangular phase diagram is mirror symmet¬ 
ric with the axes of /a = fc- Therefore one-half of the whole triangular phase diagram should 
be calculated. The triangular phase diagram obtained by our SCFT simulations is presented in 
Fig.§ Sixteen structures, including 2D, 3D ordered phases and disordered phase (D), are pre¬ 
dicted to be stable in the phase diagram. Besides the [6.6.6], [8.8.4], [8.6.6; 8.6.4], [3.3.4.3.4], 
[10.6.6;10.6.4;8.6.6;8.6.4], [12.6.4], L+C, TPL, HHC, HDG and LAM phases discussed above, 
four more ordered structures, HPL, core-shell structures of HC, DG, and BCC, and disordered 
phase are included in the triangular phase diagram. The regions of stability of the different phases 
are obtained by comparing the free energy of these candidate structures. The phase boundaries are 
determined by calculating the cross over point of the free energies of the two neighbouring phases. 
The most significant feature of the triangular phase diagram is the rich phase behavior with a large 
number of stable ordered phases. It should be noted that the ordered structures emerging in diblock 
copolymers, i.e., no core-shell cylinders, spheres, gyroid, two-color lamella, are not included in 
our calculations. That is, near the boundary of the triangular phase diagram, simulations are not 
carried out in our study. 
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a: [6.6.6] 
b: [8.8.4] 
c: [12.6.4] 
d: [3.3.4.3.4] 
e: [8.6.6; 8.6.4[ 
f: [10.6.6; 10.6.4; 8.6.6; 8.6.4] 
g: L+C 
h: HPL 
i: TPL 
k: HHC 
p: LAM 


B 



Figure 6: The triangular phase diagram of ABC star triblock copolymers with equal statistical 
segment lengths for each block and with asymmetric interaction parameters of XacN = 50.0, 
XabN = XbcN = 30.0. The ordered structures emerging in diblock copolymers, i.e., no core-shell 
cylinders, spheres, gyroid, two-color lamella, are not included in our calculations. That is, near the 
boundary of the triangular phase diagram, simulations are not carried out in our study. 
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When the chain lengths of three arms are close to each other, junction points are aligned on a 
straight line, and hence the system tends to form polygonal tilings to get smaller internal energy. 
As a consequence, the central region of the phase map is dominated by the 2D tiling patterns. The 
coordination number of each domain is proportional to the composition of corresponding block. 
For example, as discussed above, along the phase path of /a : /s : /c = 1-0: 1.0: x, the coordination 
number of C-domains changes from 4, to 6, to 8, then to 10, and finally to 12 with increasing of 
X from 0.435 to 2.76. Besides the cylindrical structures observed by experiments, two tilings of 
[8.6.6; 8.6.4] and [10.6.6;10.6.4;8.6.6;8.6.4] are predicted as the stable phases in our theoretical 
calculations. The stability has been analyzed in the above context. It is very different from the 
phase behavior of ABC linear triblock copolymer melt in which the lamellar phase occupies the 
large central region The reason is attributed to the topology of the star copolymer chain for 
which the junction points tend to be aligned on a straight line when the chain lengths of three arms 
are comparable. As the compositions of star copolymer become asymmetric, more 2D and 3D 
phases, including LAM, HC, DG, BCC, HPL, TPL and cylinders-in-lamella of L+C, appear and 
surround these 2D tilings in the triangular phase diagram. Among these phases, LAM phase has 
the following phase transitions. Near the AB edge, CACB layers can be formed whereas ABAC 
layers formed near the BC edge, and BABC layers formed near the AC edge. 

Consider the structural evolution that starts from the central 2D tilings region, toward the A- 
and C-rich comers of the triangular phase diagram. The phase transition sequence is the cylinders- 
in-lamella structures (including HPL and TPL), L-i-C, LAM, DG, HC, BCC and disordered phase. 
Near the A-rich comer, in the cylinders-in-lamellar structure, the C arms form internal cylinders in 
the B perforated lamellar, then blocks B and C form lamellar together with A-layers. The packed 
manners of the cylinders determine the morphologies of HPL and TPL, as shown in Fig.[^ The 
2D hierarchical pattern, L-i-C, where arms B and C form alternating cylinders in the lamellar-based 
phase, has a stable region between cylinders-in-lamellar structures and LAM phase. The prediction 
of stability region of L-i-C phase is generally consistent with the recent experiment in which Nunns 
et al. observed the phase in I 0 . 20 S 0 . 14 F 0.66 copolymer system Because of the smallness of 
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Xbc^ a further segregation between the minority components B and C (A and B) will not occur, 
which implies that the 3D hierarchical structures will not be global stable in this region. Instead, 
the LAM phase dominates this area. Owing to Xbc = Xab, there is similar phase behavior near the 
C-rich corner when exchanging the position of arms A and C. 

The BCC, HC, and DG phases form continuous areas across the A- and C-rich corners of the 
triangular phase diagram, where A or C block is the largest arm. The continuous areas formed by 
these phases in the triangular phase diagram reflect a continuous evolution in the compositions of 
the core and shell blocks. For example, consider the evolution of structure along a path within the 
BCC phase in the A-rich comer, starting from the AB edge, where the structure contains B spheres 
in an A matrix, to the AC edge. As the C arm increases its length, a spherical “shell” of C grows 
in the middle of each B-sphere, while the composition of the surrounding shell of B shrinks, until 
a structure of C spheres in an A matrix is obtained at the AC edge. An analogous change in the 
volume fractions of the “shell” and “core” components occurs in the HC and DG phases in both 
the A- and C-rich corners. 

The phase behavior in the B-rich corner is more complicated. Owing to the largest value of 
XacN, the minority components A and C tend to a further segregation near the B-rich corner, which 
leads to the formation of the hierarchical structures. From the cylindrical regions to B-rich comer, 
the phase sequence of hierarchical morphologies is L-i-C —)■ HDG —> HHC. Near the = fc line, 
the stable regions of HC and DG are separated by hierarchical structures of HHC, HDG, L-i-C and 
three Archimedean tilings of [12.6.4], [3.3.4.3.4], [8.8.4]. The stable region of BCC phase in the 
B-rich comer is continuous above that of HHC structure. Consider the structural evolution along a 
path that starts from the AB edge, where the length of arms satisfies the relationship /b > /a > /c> 
toward the /a = fc isopleth. Along this path, for example, in the BCC phase, A arms segregate 
into spherical domains, surrounded by C-rich pockets shell within the B matrix. The core-shell 
HC and DG phases have similar phase behavior with A-core, C-shell and B-matrix. Similarly, 
structures, such as DG, HC, and BCC, evaluating along a path that starts from BC edge toward 
the /a = fc isopleth, are of the C-core, A-shell and B-matrix patterns. In addition, the disordered 
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phase emerges in the B-rich eorner of the triangular phase diagram. 

In general, the theoretieal phase behavior is in agreement with the experimental observa¬ 
tions. However, there are some diserepancies between our theoretieal eomputationus and exper¬ 
imental observations. Some struetures observed by experiments in ISP star tribloek eopolymers, 
sueh as L-i-S phaseZine-blende type strueture^^ are not obtained in our simulations. The sta¬ 
ble regions of some struetures obtained by experiments are not exaetly loeated in the predieted 
phase regions of our phase diagram. Our theoretieal phase diagram predieted [8.6.6; 8.6.4] and 
[10.6.6;10.6.4;8.6.6;8.6.4] as stable phases which have been not observed in the experiments. 
There are three main possible reasons for these discrepancies. The first one is that in the ex¬ 
periments many star triblock copolymers are obtained by blending two kinds of copolymers, or 
adding additional homopolymers. For example, the Zinc-blende type structure was observed in 
the I1.0S2.3P0.8 system which was realized by blending the S homopolymer to Ii.oSi.gPo.s- The 
homopolymer S can definitely affect the phase stabilities. The second one is that the predicted 
phases in the theoretical calculation, such as the [8.6.6;8.6.4] and [10.6.6;10.6.4;8.6.6;8.6.4], may 
be metastable in the recent experimental systems. The third one is that there are many differences 
between theoretical systems and experimental conditions, such as the different interactions, differ¬ 
ent monomer sizes. In our calculations, we neglect the differences between XabN and XbcN. The 
difference is small, however, it might influence the self-assembling behavior. At the same time, 
the values of interaction may be different from the experimental systems. 

Conclusions 

In this work, we have investigated the phase behavior of ABC star triblock copolymers with asym¬ 
metric interaction parameters using the SCFT. Based on the previous work of screening initializa¬ 
tion strategies ,l^we can obtain a large number of ordered phases in studying the complex polymer 
systems. Then we used a fourth-order pseudospectral method, combined with the Anderson mixing 
algorithm to calculate the free energy of the observed phases. Motivated by previously experimen- 
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tal studies, the Flory-Huggins interaetion parameters of Xbc^ = Xab^ = 30.0, Xac^ = 50.0 are 
used to model the experimental systems of ISP and ISF star eopolymers in bulk. To extend the 
seope of theoretieal study, a large number of 2D and 3D ordered struetures have been involved in 
our ealeulations. In order to shed light on the phase behavior of the ABC star tribloek eopolymers, 
we first determined the phase stability along the phase path /a'- /b '■ fc = 1-0 : 1.0 : .r. Our results 
agree with those of ISP star triblock copolymers for cylindrical phases well. Then we calculated the 
phase regions along the phase path /^ : /b : /c = 1.0 : .r: 1.0. On this phase path, besides the cylin¬ 
drical structures, we emphatically analyzed the stability of hierarchical structures of L-i-C, HDG, 
HHC. The phase stability has been analyzed by splitting the SCFT energy functional into internal 
and entropic parts in detail. Finally we constructed a very complicated triangular phase diagram 
with these candidate structures. Owing to the case of interactions XbcN = XabN, the phase diagram 
has only one mirror symmetric axis fA = fc- Fifteen ordered structures and the inhomogeneous 
phase constitute the triangular phase diagram. In general, the phase regions predicted by our SCFT 
calculations are consistent with previous theoretical studies of either SCFT calculations or MC, 
DPD simulations which mainly focus on the equal interaction systems. However, in our calcula¬ 
tions, the interaction parameters are more closely related to experimental systems, such as ISP and 
ISF star triblock copolymers , i.e., the asymmetric interaction parameters XacN > XbcN = XabN. 
It has been found that the asymmetry of the interaction parameters plays a profound role in the 
complex phase formation. Furthermore, our predicted phase diagram involves more phases, es¬ 
pecially 3D structures, and presents more comprehensive phase behavior. Our calculations and 
analysis can be helpful to understand the self-assembling mechanism of complex structures. The 
resulting phase behavior extends the theoretical study to the asymmetrically interacting ABC star 
triblock copolymers. And the presented phase diagram will be a useful guide for further study of 
ABC star triblock copolymers. 
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Appendix: The Implementation of Anderson Mixing Method 


The k iteration in the Anderson mixing method begins with the evaluation of new fields from 


Eqns. (??)-(??) 


(fc). ^ ^ La=A,B,cw^a (r)^a " ^NxabXbcXac 
“ LaXa 

w^A (*■) = XabN^P (r) + XacN^^c^ (r) + r](r), 
w^B (r) = XABN^f^ (r) + XbcN^P (r) + t] (r), 
wf (r) = XAcN^f^ (r) + XbcN^^b^ (r) + t] (r). 


( 12 ) 

(13) 

(14) 

(15) 


In the above expressions, a =A,B,C, are the old fields, Xa = Xbc{Xab + Xac-Xbc), Xb = 
XAciXBC + Xab - Xac), Xc = Xab(Xac + Xbc - Xab)- Next we evaluate the deviation. 


jW = w^^((r) — w(^^(r), 


(16) 


where = {df^ {r),dP {r),dP {r)Y, w(*((r) = (w^^^(r), w5fV)5M;?(r))^, w(*((r) = {wf\r), 
(r), ^ (r) ) ^ • From the deviation we can specify an error tolerance through the inner product 


(^(r),/(r)) = ^ j drg{r)f{r), 


(17) 


where g(r) and /(r) are arbitrary functions. The error tolerance is defined by 


ei 


' Za=A,BMa\r)A\r)) 


1/2 


(18) 


as a measure of the numerical inaccuracy in the field Eqns. (??)-(??). 

The simple mixing method is performed until a certain tolerance is reached where a morphol¬ 
ogy has begun to develop. From our experience, Ci < 10^^, is sufficient in most cases. We then 
switch to the Anderson mixing procedure by the previous n steps to update fields. We assemble 
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the symmetric matrix in this way 




(19) 


for ij— 1,2,..., n, and vector 

y{k) ^ j(*)) _ (20) 


From these, we calculate the coefficients 


Ci = E(t/-‘)ijV,-, (21) 

j 


and combine the previous histories as 


I(x — W(x + 

!=1 

(22) 

i=\ 

(23) 


Finally, the new fields are obtained from 


(A:+l) 


r« + AD 


(k) 

ccj^ 


(24) 


where 0 < A = 1.0 — 0.9^ < 1 

In our implementation, the used previous steps are usually much less than the number of basis 
functions or grid points. Therefore assembling the n-order linear system spends more computation 
time than solving this system. To save computational amount, we decompose Eqn. (??) into 

ulf = (jW, jW) - jW) - (25) 

In k iteration, only the terms related to in the right term of Eqn. (??) are required to calculate. 
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but the last terms i,j= 1,... ,n, should be not computed repeatedly which will 

save the main computational cost in assembling matrix U. In practical calculations, we store the 
following inner product matrix 




0 


( 26 ) 


Then we can assemble the matrix U and vector V according to the expressions of Eqns. (??) and 
(??) using the elements of matrix S. Note that the inner product matrix is symmetric, therefore 
only a row (or a column, equivalently) of 5 is required to update in each iteration. 
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